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Introduction Générale : 


L'utilisation des explosifs en génie civil a d'importantes applications 
depuis des siècles, e.x. : extraction de roche dans la construction de 
voies de communication, excavation pour tunnel souterrain, démolition 
de bâtiments, exploitation minière … 


La conception de ces opérations, dites tirs de mine ("rock blasting") est 
une tâche délicate. Il s'agit de choisir le type d'explosif, sa distribution 
géométrique et les différentes façons de déclencher les tirs pour une 
géologie donnée. Le meilleur résultat des tirs pour un volume de roches 
abattues a pour objectif la qualité de la fragmentation pour un coût de 
revient donné sans violer les conditions de sécurité ; c'est là, une 
solution optimisée. En raison de la complexité du problème, les 
techniques des tirs sont restées durant de nombreuses années un art où 
l'expérience du praticien prend une part dominante. Les ouvrages 
classiques [LK63],[Hin59], [Hem]... montrent des expériences dans 
lesquelles les formules empiriques sont les principaux outils à la 
conception. Les facteurs qui influencent le résultat des tirs sont 
nombreux. Nous citons ci-dessous les principaux : 


e caractéristiques des explosifs pour le dégagement d'énergie 
explosive 

e caractéristiques mécaniques des roches 

e caractéristiques du plan de tirs qui décrivent la géométrie de 
foration 

e caractéristiques du système d'initiation des tirs. 


L'expérience montre qu'une distribution des explosifs dans le massif 
rocheux avec les temps de retards d'initiation conduit à un meilleur 
résultat 


Du point de vue mécanique, l'utilisation des explosifs dans l'excavation 
d'un massif rocheux implique deux branches théoriquement difficiles : 
la première sur la mécanique des explosifs, i.e. la détonique ; et la 
deuxième sur la mécanique de roches en dynamique. Leur difficulté 
commune est le couplage des problèmes locaux/globaux. Pour la 
première, il s'agit du couplage physique/chimique fortement non 
linéaire; et pour la deuxième, il s'agit de tenir compte des défauts 
(notamment des fissures) du milieu rocheux (la modélisation du 
comportement statique pour un simple essai en compression est déjà 
assez compliquée). L'étude sur ces deux branches est encore en cours 
de développement. Néanmoins, la prévision fondée sur les modèles 
mathématiques des résultats des tirs est une nécessité en raison de 
motifs économiques évidents (les essais expérimentaux ainsi que 
l'acquisition des informations sont coûteux et les modèles quasi- 
empiriques ne sont valables que pour des problèmes spécifiques). Les 
résultats analytiques de problèmes de cavités [Sha42], [Duv50], 
[Fav69]... sous certaines hypothèses bien fortes étaient la principale 
base de la prévision rationnelle. Depuis ces dernières années, des 
progrès considérables, liés à la puissance des ordinateurs ainsi qu'aux 
traitements numériques avancés permettent de développer des outils 
(simulateurs) plus sophistiqués pour la prévision de la qualité, ou bien 
encore pour approcher une conception optimisée. 


JUSTIFICATION DU CHOIX : 


L'utilisation des explosifs en génie civil a d'importantes applications 
depuis des siècles, e.x. : extraction de roche dans la construction de 
voies de communication), excavation pour tunnel souterrain, 
démolition de bâtiments, exploitation minière … 


PROBLEMATIQUE : 


Du point de vue mécanique, l'utilisation des explosifs dans l'excavation 
d'un massif rocheux implique deux branches théoriquement difficiles : 
la première sur la mécanique des explosifs, i.e. la détonique ; et la 
deuxième sur la mécanique de roches en dynamique. 


Donc comment peut-on remédier mettre fin à ce problème ? 


SOLUTION : 


application des progrès considérables, liés à la puissance des 
ordinateurs ainsi qu'aux traitements numériques avancés permettent 
de développer des outils (simulateurs) plus sophistiqués pour la 
prévision de la qualité, ou bien encore pour approcher une conception 
optimisée. 


2. L'optimisation : 


L'optimisation concerne la base même de la conception. Dû à une 
grande difficulté d'utilisation, les méthodes d'optimisation sont très 
variées, chacune conçue pour résoudre une catégorie particulière de 
problèmes. Elles sont développées souvent indépendamment les unes 
des autres. Par exemple, l'évolution des méthodes d'optimisation des 
structures a commencé par l'intuition (e.x. la tour Eiffel est vérifiée 
comme étant une structure optimisée), suivi par des règles simples (e.x. 
"Fully Stressed Design"). Depuis une vingtaine d'années, liée au succès 
de l'utilisation de l'ordinateur dans le domaine de l'analyse des 
structures discrétisées (éléments finis, éléments frontières...) la 
méthode dite Séquentiel Quadratique Programmation (SQP) est le plus 
souvent utilisée dans ce domaine. 


On classifie ces divers techniques en deux grandes catégories; la 
première concerne des méthodes traditionnelles (dont SQP fait partie) 
et la deuxième concerne les méthodes d'évolution. 


2.1 Méthodes traditionnelles : 


Guidées par un concept de continuité, les méthodes traditionnelles sont 
originaires des applications de la programmation mathématique 
(linéaire, non-linéaire) fondée sur la résolution des conditions 
d'optimalité de Kuhn-Tucker, pour lesquelles les méthodes de type 
gradient sont essentielles. 


La formule mathématique générale pour une fonction objectif / s'écrit : 


trouver min f(x) où B=f{x:gx) 0; i=1,... Ng} 


qui consiste à chercher un vecteur des variables (dites variables de conception) x € R", 
tel que la fonction de coût (fonction objectif) à minimiser (ou à maximiser) f : R" => R 
sur un domaine admissible B défini par les fonctions de contrainte g; : R" = R 

Pour une programmation non linéaire, les procédures de résolution sont fondées en 
général sur les 1térations 

F+! = 2 4 opt où p'est une direction calculée à partir des gradients de la fonction 
objectif et des contraintes actuellement satisfaites (ensemble actif à l'itération $) ; a, pas 
de descente doit être évalué dans une autre étape. 


De nombreuses techniques appelées méthodes de descente (plus 
grande pente, gradients conjugués, quasi-Newton, recuit simulé etc..) 
ont été développées afin d'augmenter l'efficacité de l'algorithme ; pour 
prendre en compte que l'espace de recherche est de très grande 
dimension, que la fonction de coût et les contraintes ne sont pas 
convexes et que les contraintes sont impossibles à exprimer 
quantitativement. Ces techniques n'apparaissent malheureusement que 
sur certains problèmes et pas sur d'autres, leurs faiblesses sont 
généralement sur : 


e Minima locaux 
la solution est piégée facilement dans des minima locaux 
e Hypothèses restrictives 
les fonctions f, g sont continues et dérivables 
e Le traitement des différents types de données 
elles ne traitent que des variables continues 
e Spécificité 
Les techniques d'optimisation ne sont jamais générales. 
2.2 Méthodes évolutionnistes : 


En rupture avec les méthodes traditionnelles, pour balayer l'espace, les 
méthodes évolutionnistes font "évoluer" une population de solutions 
d'une façon suggérée par la génétique. 


On donne les quatre concepts particuliers qui distinguent les méthodes 
évolutionnistes des autres méthodes d'optimisation : 


1. Population Une méthode évolutionniste entretient une population 
d'individus ("population de solutions"), chaque individu est 
comme un point du vecteur des variables de conception dans la 
méthode de descente. Il ne s'agit pas simplement d'effectuer 
plusieurs recherches en parallèle à partir de plusieurs solutions 
différentes, le concept de population est essentiel, car les 
différentes solutions de la population de la solution vont interagir. 

2. Evolution 
La nouvelle population est créée : 

e soit par modification d'un individu ("mutation") 

e soit par interaction entre plusieurs individus déjà existants 


On dira que ces nouvelles solutions (individus) sont les "enfants" 
des solutions qui ont interagi, et que celles-ci étaient leurs 
"parents". 

3. Valeur ("fitness") Cette valeur pourra, dans les cas les plus 
simples, être simplement le coût associé à une solution (une 
valeur de la fonction objectif dans la méthode de descente), mais 
pourra également être une grandeur plus complexe ( la valeur 
d'une fonction généralisée qui contient les termes de pénalité des 
fonctions de contraintes). 

4. Sélection Pour garder une population constante après ou avant 
l'évolution, la sélection consiste d'éliminer les individus considérés 
moins bons que les autres. Cette élimination est soit probabiliste, 
soit déterministe selon les particularités de la méthode 
considérée, mais elle est essentielle aux méthodes 
évolutionnistes. 


3. Optimisation couplée à l'apprentissage 
automatique : 


L'optimisation des tirs à l'aide des méthodes traditionnelles de 
l'optimisation apparaît impossible. D'abord, dans la conception des 
tirs, il y a beaucoup de variables de conception qui ne sont pas du 
type continu (e.x. modes d'amorçage des explosifs), ensuite, les 
simulations sont trop coûteuses pour suivre les itérations. En 
revanche, les méthodes évolutionnistes comme les algorithmes 
génétiques sont bien adaptées pour résoudre des problèmes 
d'optimisation généraux, mais il faut compenser leurs faiblesses 
notamment le nombre d'itérations et le temps des calculs. 
L'apprentissage automatique permet de satisfaire ces besoins, car il 
peut permettre de générer des règles (fonctions) nécessaires qui 
seront intégrées ensuite dans les algorithmes génétiques. Il faut 


quand même résoudre un problème très délicat, la construction de la 
base des données soit par simulation numérique, soit par essais, soit 
encore en couplant des résultats analytiques, numériques et 
expérimentaux. Nous n'avons pas eu la possibilité d'effecuer des 
essais avec explosifs durant cette thèse. Ce chapitre montre 
cependant comment effectuer la génération des règles par 
apprentissage automatique et le couplage à l'optimisation à partir 
d'une base de données construite par simulation numérique. Les 
exemples traités sont certainement loin de la réalité, mais servent de 
démonstrations. 


3.2 La recherche des règles par l'apprentissage 
automatique : 


L'option de normalisation du type X = X / max [abs X) est choisie pour 
tous les exemples effectués 


3.2.1 Coefficient polytropique : 


Il s'agit de chercher une expression simple du coefficient 
polytropique 7 en fonction des caractéristiques des explosifs 
notamment a, le rapport énergie de gaz/énergie de choc obtenu lors 
de l'essai en immersion. D'après les conséquences (détonation dans 
l'eau), il y a trois caractéristiques indépendantes dans une loi 
polytropique, les autres peuvent être 
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Tableau 15.1 La base de données du coefficient polytropique 
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estimées par les formules simplifiées. On a choisi les quatre descripteurs suivants : 

densité initiale d’explosif p. , énergie dégagée E, 7. et a. 

La base de données est construite par des simulations numériques similaires à celles 
indiquées dans le chapitre 5. 

16 calculs ont été effectués et ont conduit au tableau (15.1). 

Les valeurs maximales des descripteurs sont : 

max(p.) = 1.0;max() = 3.0; max( E) = 0.01155 ; max(a) = 2.313. 

En prenant + comme conclusion, SEA donne le résultat suivant : 


è avec les options : polynôme d'ordre 1 et les termes rectangles, 


4 = 1.58 + 4.0195E, — 2.256E? — 2.045E, : a (15.1) 


(la moyenne de la somme du carré des écarts sur le test est €m = 0.016.) 


e avec les options : polynôme d'ordre 1, les termes rectangles, et les termes inverses : 


1 1 
y = 0.466—— — 0.1 -=z + 1.54 15.2 


(et €m = 0.00106.) 
On prend ensuite œ comme la conclusion. SEA donne le résultat suivant : 
e avec les options : polynôme d’ordre 1, les termes rectangles, et les termes inverses : 


a= 2.74% — 1.9176p. (15.3) 


(et €m = 0.0368.) 


e avec les options : polynôme d'ordre 2, les termes rectangles, et les termes inverses : 


E 
a = 1.96 - 1.88467 + 0.9977 + 0.978ÿE? — 1.55454 (15.4) 
| 


(et em = 0.001522.) 


Dans les apprentissages ci-dessus, on remarque que Er, pe ne sont jamais présents en 
même temps ; cela est dû au fait qu'il existe une relation de dépendance entre ces deux 
descripteurs dans la base de données. SEA a pu éliminer un de ces deux descripteurs 
dépendants. Les formules (15.2),(15.3) semblent plus simples et plus logiques (physique- 
ment). La formule (15.3) est enfin la plus pratique, puisque E; est obtenue souvent par 
estimation. 


3.2.2 Loi de propagation : 


Une loi classique pour prédire l'intensité des ébranlements par les tirs à explosif, dite loi 
de propagation est la suivante : 

KË 

ra 


v=K 


où 

v est la vitesse particulaire résultante, Q est le poids d’explosif utilisé, d est la distance 
au tir et K, a, b sont des constantes caractérisant l'énergie explosive et la roche. 

En effet, cette formule très simple peut être déduite par analyse dimensionnelle : tableau 
(9.1) de la section 9.2, le rayon de la cavité est proportionnel à Qt pour la cavité sphérique 
et Q? pour la cavité cylindrique. 

Comme cela avait pu être noté dans le chapitre 12, l’accélération est une caratéristique 
du niveau de vibration commode grâce sa propriété d'indépendance de la rigidité et du 
rayon de cavité. On a choisit donc les descripteurs suivants : 





conclusion pr a(—#—) | produit de l'accélération par la densité du rocher | 
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Pour construire la base de données, 9 calculs élastiques en axisymétrie ont été effectués 
sur la géométrie donnée dans l'exemple 3 du chapitre5; pour chaque calcul les accélérations 
aux 6 points différents sont enregistrées. 

Les valeurs maximales absolues des descripteurs sont : 

max( E) = 73400 ; max(Ÿ) = 3.0 ; max(E;) = 1550; max(v) = 0.401 ; max(d) = 814.9. 

Avec les options : polynôme d'ordre 1, et les termes rectangles, SEA donne la conclusion 
suivante : 


pra = fo — fa d + 4.57d? (15.5) 


où : fo = —2.39 + 3.45E, + 1.29% - E — 0.19v* + 5.22v : Ÿ + 1.597 - E, + 0.78E* 
fa = —3.88 + 1.21E + 4.53v + 6.347 + 5.88E, 

(et en = 0.00097.) 

Si toutes les valeurs des descripteurs sont prises en logarithme, les valeurs maximales 
absolues des descripteurs actuels sont : 

max(E) = 11.204; max(Ÿ) = 1.0986 ; max(E;) = 7.34; max(v) = 1.772 ; max(d) = 6.25. 

En choisissant les options : polynôme d'ordre 1 sauf pour v qui est d’ordre 2, la con- 
clusion approchée est la suivante : 

Log(p,a) = —9.976 + 56.8 Log v + 30.72(Logv)? + 30.43 Log E + 19.68 Log E; + 4.38 
Log Ẹ — 16.82 Log d 

(et em = 0.01.) 

Le module d'Young E intervient dans l'apprentissage, c'est parce que la pression de 
l'explosion dépend de cette grandeur (voir équation (9.18) de la section 9.2) dans le calcul 
du couplage détonation-structure. 


3.3 Exemple d'optimisation : 


Le problème posé est le suivant : 

Pour un type d’explosif donné, un type de roche donné, on cherche les valeurs des 
paramètres du plan des tirs (d,R;,B) qui maximisent l'efficacité spécifique 

E, = 


2 
Nous reprenons le résultat de SEA (l'équation (15.7) ) divisée par (4) qui est pro- 
portionnelle au volume explosif. Le problème s'écrit mathématiquement : 


RER RE 
ve(volume de l'explosi£ utilisé) ' 


2 
max f = 12568% — 88.95B + (90.39 — 92.644R4 + 25 HR) EE 


avec les contraintes : 

1.0 < Ri < 2.0 

04< B<1.0 

0.5 < d < 0.2 

d/ R4 > 0.03 (le diamètre de la cartouche explosive est supposée supérieur à 0.03m) 
Ce problème ”simple” est résolu d'abord par une des méthodes traditionnelles : 


e Gradient Conjugué 


Le point de départ est {d = 0.04, B = 0.8, Rd = 1.0} , une solution a été trouvée telle 
que : 
{d = 0.04, B = 0.5, Rd = 1.09358} = f = 1187.1658 


En modifiant le point de départ en {d = 0.1, B = 0.8, Rd = 2.0} june solution différente 
de la précédente a été trouvée telle que : 


{d = 0.06, B = 0.5, Rd = 1.09358} => f = 995.79 


Cela signifie qu'il existe au moins deux points de solution locale et que la solution 
trouvée dépend du point de départ. 


Ensuite, on utilise l'algorithme génétique pour résoudre le même problème : 


e Algorithme Génétique 
A partir de plusieurs points différents, une même solution a été trouvée telle que : 


{d = 0.04, B = 0.5, Rd = 1.094118} = f = 1187.1654 


qui est pratiquement la même que la première solution trouvée par la méthode de 
Gradient Conjugué. 


Maintenant, supposons que le diamètre du trou et le diamètre de la cartouche explosive 
ne peuvent prendre que certaines valeurs fixes telles que (cela est peut-être dû à la limite 
de l'équipement de foration ainsi que des explosifs disponibles) : 


d € {0.048,0.05, 0.066, 0.084, 0.1, 0.15, 0.2} 
# € {0.03,0.036,0.04,0.05, 0.64} 
Les méthodes traditionnelles ne peuvent pas traiter directement des variables discrètes 
ou multivalués ; les algorithmes génétiques sont mieux adaptés pour ce genre de problèmes. 
Voici la solution trouvée par l'algorithme génétique en changeant la variable Ra par 
d/ Ri : 


{d = 0.048, B = 0.5, d/Ra = 0.045 > Rd = 1.0667} = f = 980.96 


R.q: Le calcul avec une méthode du type gradient est naturellement plus rapide qu'avec 
l'algorithme génétique, c’est la raison principale pour laquelle les méthodes tradi- 
tionnelles sont encore dominantes dans le domaine de l'optimisation des structures 
où le nombre de variables de conception est souvent important. Mais, l'algorithme 
génétique couplé avec l'apprentissage automatique permet de compenser l'inconvénient 
tout en gardant ses avantages. 


3.2.3 Volume endommagé : 


Il s’agit cette fois pour un explosif donné et une roche donnée (fixes), d'approcher une 
expression du volume endommagé va(roche abattue) en fonction des caractéristiques du 
plan de tirs afin d'effectuer le calcul d'optimisation du plan de tirs. 

La base de données est construite par des calculs en déformation plane dont les config- 
urations sont similaires à celle dans l'étude de l'influence des macro-discontinuités (figure 
13.1 ). Il y a trois surfaces libres, si l’on compte aussi la symétrie de la configuration. Les 
descripteurs choisis sont : 

diamètres du trou d ; distance entre le trou et la surface horizontale B; (la distance 
entre le trou et la surface verticale est fixée à 0.75), rapport de découplage R4 ainsi que 
deux rapports d/R, et d/B. 

Onze calculs ont été effectués. Cela donne les résultats dans le tableau (15.2). 


iom a B STEEL EE 
eame 00s0) 0.600 | 10 [192774] 005 100833] 
l ed1bac3 0.050 0.5066 | 1.677 | 0.79338 | 0020 | 0.0087 | 
Cealbrez | 0.050 | 0.600 | 1.172 | 1467875 | 0.0427 [0.0833 | 
[edībic3 | 0.050] 0.600 | 1.667 | 0.673639 | 0.0298 10.0833 | 
[edzbic1 | 0040| 0.500 11.3333 | 094768 | 003 | 008 | 
“ed2bic2 | 0.010 | 0.000 | 1.333 [0871909] 0.05 10.0667 | 
Ced3b2e2 | 0.100 | 0.750 | 1.257 | 1.320870 | 00798 | 0.133 | 
Led3b2c3 [0.100 | 0:750 | 1.474 [1266204 | 0067 | 0.133 | 
Code? [0.200 | 1000 | 1.2067 | 1570654101657] 02 | 
[edab2e2 | 0.200 | 0.750 | 1.2067 | 2.957848 | 0.1657 | 0.267 | 
Codes [0.200 | 1.000 | 1.843 | 1182392 [0108 | 02 | 


Tableau 15.2 La base de données du volume endommagé 






















Les valeurs maximales des descripteurs sont : 
max(d) = 0.2;max(B) = 1.0; max(R4) = 1.843; max(d/Ry) = 0.1657; max(d/B) = 
0.267. 


Avec les options : polynôme d'ordre 1 et les termes rectangles, SEA donne un résultat 


simple : 


d d 
= 8.1— -0.65B— 
Vd Pi 65 Ru 


(et es = 0.1.) 


Avec les options : polynôme d'ordre 2 et les termes rectangles, SEA donne le résultat 


suivant: 


2 
va = 0.85B : d — 2.44B ($) + (24.134 — 45.6R4 + BTRÈ)E 
d 


(et em = 0.00477.) 


Cette formule est retenue dans l'étude d'optimisation qu'on montrera dans la section 


suivante. 


Application au Maroc : 


Dans le Maroc, comme dans le monde entier, L'utilisation des explosifs 
dans les travaux de génie civil est très diversifiée : elle est destinée aux 
terrassements rocheux de voies de communication à ciel ouvert ou en 
souterrain (routes, voies fluviales, voies ferrées), aux creusements de 
fondations d'ouvrages d'art (barrages, ponts et centrales) et de 
bâtiments, aux aménagements portuaires, à la démolition et à certains 
emplois spéciaux, notamment dans la sidérurgie. 
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